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Abstract 

We consider the classical equations of motion for a single Galileon field with generic 
parameters in the presence of non-relativistic sources. We introduce the concept of 
absolute stability of a theory: if one can show that a field at a single point — like infinity 
for instance — in spacetime is stable, then stability of the field over the rest of spacetime 
is guaranteed for any positive energy source configuration. The Dvali-Gabadadze- 
Porrati (DGP) model is stable in this manner, and previous studies of spherically 
symmetric solutions suggest that certain classes of the single field Galileon (of which 
the DGP model is a subclass) may have this property as well. We find, however, that 
when general solutions are considered this is not the case. In fact, when considering 
generic solutions there are no choices of free parameters in the Galileon theory that will 
lead to absolute stability except the DGP choice. Our analysis indicates that the DGP 
model is an exceptional choice among the large class of possible single field Galileon 
theories. This implies that if general solutions (non-spherically symmetric) exist they 
may be unstable. Given astrophysical motivation for the Galileon, further investigation 
into these unstable solutions may prove fruitful. 



1 Introduction 
1.1 The Galileon 

The Galileon naturally arises when one considers ghost free theories that modify General 
Relativity in the IR [T] (see also, [2]-[Z! and its covariant extension [8J and [9J, and its 
use as an alternative to inflation [ID])- In a local patch smaller than the cosmological 
horizon the additional degrees of freedom that encode our new dynamics are given by a 
relativistic scalar field, denoted by vr, which is universally coupled to matter via vrT^, where 
71 generically decouples from matter at short scales due to derivative self-interactions (dubbed 
the Vainshtein effect). The Lagrangian describing these extra degrees of freedom is given by 

5 

Ct, = M^iCjCj , where schematically £„ = {d'^7!-)'^~'^dndn (1) 

i=l 

and the equation of motion for the field is 

E^^^^ = -17^' ^^^^^ £n = ^ = {dd7rr-' (2) 
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We arrive at this particular equation of motion by imposing the symmetry 71 (x) — ?■ 
7r(x) + b^x^ + c. This symmetry will require the equations of motion to have at least two 
derivatives. However, a well defined Cauchy problem/absence of ghosts requires that the 
equation of motion have at most two derivatives. These requirements are strong enough to 
uniquely determine each (these can be found in [1]) to an overall constant, given here by 
the Cj's. This structure is very powerful, as we can re-write the partial differential equation 
as algebraic one 

F{ki_,i,) = -p , where kf,„ = d^dy-n and p = Tj^ (3) 

In order for the Galileon to be an interesting modification to gravity (its original motiva- 
tion), we demand that the non-linear nature of the Galileon theory becomes manifest at the 
cosmological scale; that is, we want the non-linearities to become important at Hubble scales. 
This allows, in principle, cosmological solutions driven by vr's self-interactions rather than by 
the usual matter T^,^. We rewrite CiCi ~ Cidnd-w {O'^ti/Hq)"' where the new dimensionless 
coefficients Cj, defined as the combination of original ones Cj and Hubble constant Hq, should 
be of order unity. A quick word of clarification: when we make reference to the DGP model, 
we mean not strictly the coefficients that come from the 4-d effective field theory associated 
with the extra-dimensional model put forth by [H], but rather the class of theories where 
ci = 0, £2 > 0, £4 = 0, and £5 = 0. Subsequently, Galileon terms have also appeared in 
the decoupling limit of an interesting interacting theory of Lorentz invariant massive gravity 
[13] as well as in a supersymmetric extension of ghost condensate theories \1A\. Additionally, 
the following analysis could possibly be extended for the more general cases of multi-field 
Galileo ns P^ . 

In this paper, we are concerned with classical stability and the non-linearities perform 
a critical role. We can investigate such a classical question because we are assured that 
quantum effects are small. Thanks to the Galileon invariance, the £nS above are not renor- 
malized under loop corrections, and terms with fewer derivatives acting on the vr field are 
not generated quantum mechanically [12] (a similar proof holds for the multi-field Galileons 
[15]). In particular, upon estimating the one loop effective action one can see that for clas- 
sical solutions with large non-linearities quantum effects are small provided gradients are 
mild 9 <^ A [16j. All of the above indicates that it is consistent to focus on the stability of 
classical solutions to 

In [T7] it was shown that there do not exist any stable, static, source free soliton solutions. 
Any such solution, if it existed, would be extremely unstable. This instability lies in 'wrong 
sign' kinetic terms, and therefore instabilities exist all the way down to the UV cutoff. 
Such a solution, if you attempted to construct it, would simply exit your effective theory 
immediately. This should be contrasted to instabilities stemming from 'wrong sign' potential 
terms. In this case one can follow a constructed solution as it changes into something else 
while still maintaining perturbative control (like a ball slowly rolling down shallow hill). 
We ask the following: does the picture change for non-linear solutions when non-relativistic 
sources are added to self interactions? Are these classically stable? Some analysis has been 
done along these lines. In particular, radial solutions were throughly investigated in [1]. 
Given the recent interest in Galileon theories, we feel it prudent to build upon, generalize, 
and possibly clarify the work done for radial solutions. 
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This paper is organized as follows: In Section 2, we specify our problem and describe 
the general logic of our analysis, followed by a brief report of our results. In Section 4, in 
order to make our results more transparent, some graphical and analytical interpretations are 
presented. Detailed (and lengthy) analysis is contained within Section 3 and the Appendix. 
Readers are invited to skip these sections as they have been designed to contain all the gory 
details. In Section 5, we summarize our findings and suggest avenues for future inquiry. 

2 Outline of Analysis 

We begin by assuming some solution, say ttq, for the Galileon field that satisfies the equa- 
tions of motion in the presence of some source, and then expand around this solution by 
adding small fluctuations. Taking ttq — ?■ ttq + and keeping second order in terms in the 
Lagrangian, we have the dynamics of the fluctuations!^ The action is given by 

^^ = lj d'xZ'^''{x)d,^d,^ (4) 

where Z^j, is a matrix made out of the second derivatives of our ttq field (i.e. made out of 
the matrix fc/ijy)]^ 

On the other hand, one can show that the matrix Z^^^ can be found by taking the first 
order variation of the equation of motion 

6£n = -Z^^d^dA (5) 

which provides an easier way to compute Zj^ Simply put, we will investigate under what 
conditions these small fluctuations remain small. 

2.1 Conditions for Stability 

The main goal of this paper is to explore the parameter space spanned by the q's [i = 
1 . . . , 5) to determine if there exists a subspace in which the Galileon theory is stable. We 
should point out that the "stability" we refer to is consistent with the physical meaning meant 
by [LQ}. While there are other definitions of stability, we utilize this definition because, 
though it may be limited, it is precise. Given that it is slightly different than the usual 
definition of stability in ODEs or PDEs, it is worth a few explanations. First, the stability 
we consider here is a local one, i.e., on a space and time scale much shorter than those 
typical of the background field ttq, and thus we are safe to treat the matrix function Z^jy(x) 
as constant. Therefore (it should be noted) we will not be able to keep track of phenomena 
like resonances which can be interpreted as "instabilities" . These "instabilities" are of a much 

^Throughout this paper we wiU be using the (— , +, +, +) signature. 

•^The reason why only terms proportional to d4>d(j) survive in ^ is as follows: Note that the nth order 
Lagrangian can be written as £„ = r''i'^i ' '^"'^"(9^j7r9^j7r3^t2^2 7r . . . d^^^^ir, where T is antisymmetric under 
changes of any [jii^jij) or (vi^Vj) pair while symmetric under that of {pi,Vi) pair jlj. Thus, thanks to the 
properties T possesses, all the terms in the second variation of £„ can be expressed as 909(/)(9i97ro)"~^ + 
surface terms 

•^Note that this is consistent with Q because d^Z^" = identically. 
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less catastrophic nature than the ones we discuss and, as previously mentioned, still allow 
for much analytic control. In the neighborhood of a given point, our stability corresponds 
to demanding the e.o.m. of the fluctuation field 



(6) 



give an oscillating solution. When cast in Fourier space, oscillating solutions correspond to 



having real solutions for uj for any real spatial momentum vector g, or equivalently, the 
matrix ZqiZoj — Zo^Zij will be positive definite. Second, we want the Galileon theory to be 
absent of ghost-like instability, i.e. the sign of the kinetic term of the fluctuation action ^ 
to be correct, which requires Zqq > 0. 

As previously mentioned, we will focus on non-relativistic matter sources, not only be- 
cause of their great importance, but because it can be shown that, given such sources at a 
generic spacetime point, the symmetric tensor /c^^,, and therefore Z^^, can be diagonalized 
through a Lorentz transformation [16] . Thus, the conditions for the local stability are simply 



where Z^'s are the diagonal elements of the matrix Z^ diag(Zo ,^1,^2, Z3). 

It has been shown in the DGP model that, given positive energy density sources, if 
a specific solution is stable at some point (in the way indicated above) then its stability 
throughout the spacetime is assured [16]. This is, of course, a desirable property of a theory. 
One may wonder whether a generic Galileon theory shares this same nice property. Are there 
"safe" choices of the Cj parameters (evidently, they should include the DGP parametrization) 
such that this subclass of Galileon theories possess the same property as the DGP model? 
The answer to this question is the main result of this paper. 

In order to be absolutely clear about what we accomplish, we define the concept of 
absolute stability. Assume some solution to the equations of motion, ttq, exists. 

Definition 1. An absolutely stable region in parameter space is a region of {ci, C5} 's 
where, if at a single point in spacetime, say Xq, Z^(7ro(xQ)) < (i.e. ttq is stable at this 
point), then for non-relativistic source profiles satisfying p G [0, 00) the equations of motion 
guarantee that Z^ < over the rest of spacetime (i.e. ttq is stable over all of space). 

Why is this a useful concept? If a choice of parameters is absolutely stable, then it follows 
that for any non-relativistic positive energy source configuration — no matter what the global 
structure — stability of a particular solution at a single point implies global stability. When 
talking about absolute stability one does not have to solve (or at least characterize the 
solutions) the Cauchy problem for all possible source configurations. We consider the whole 
equation of motion surface rather than characterize particular solutions. This stronger cut 

^When we formulate our analysis in terms of the eigenvalues of the k^^ tensor, the equation of motion 
£ — p/Mpj (at a single point in spacetime, or equivalently, with p fixed) defines a surface in the space spanned 
by these eigenvalues. It is this surface that we are referring to. All the surfaces generated by different sources 
{p G [0, 00) ) we group into a "family" . 
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on acceptable parameters allows us to side-step the difficulties of dealing with the arbitrarily 
complicated global structure of our source. Considering that we are dealing with a non-linear 
PDE it is surprising that we can say anything at all. In the general case we don't know how 
to show existence, but we are still able to say something about stability. 

The DGP theory is absolutely stable [16]. Spherically symmetric solutions (and mild 
deformations of them) for particular choices of parameters in a single field Galileon model 
PP are stable. Is there some part of this parameter space that admits absolute stability? 

To our surprise, we find that the DGP model is the single absolutely stable class. That 
is, the powerful property of absolute stability that the classical DGP theory possesses does 
not carry over into the general Galileon theory. 

2.2 General Program 

After diagonalizing fc^j^ with an appropriate boost we can write the stability conditions in a 
nice algebraic way: 



Zo = Zo{c2,C3,C4,C5,ki,k2,k3) <0 (9) 

Z^ I = Zi(c2,C3,C4,A;2,fc3) < (10) 

Z'^2 = ^2(C2,C3,C4,A;1, fcg) < (11) 

Z% = Zsic2,C3,Ci,ki,k2) <0 (12) 



where the q's are the coefficients that describe our freedom in choosing the exact Galileon 
Lagrangain ([T]) and the fc's are the eigenvalues of /c^^ for non-relativistic sources. The reason 
ko (= koo after the matrix has been diagnalized) does not appear in the expression above 
is because it is suppressed by two powers of f <^ 1 in comparison to the other eigenvalues, 
that is fcoo ~ v'^kij, which must be small by assumption in order to ensure diagonalization. 
Zq is a cubic function of the fc's while the Zj's are quadratic. 

Additionally, we have the equation of motion for ttq which becomes an algebraic equation 
for the fc's. Note that Ci and p enter in the same manner and are easily combined when we 
consider a single point in spacetime. We have 

S{ci, C2, C3, C4, C5, ki, ^2, ks) = (13) 

S is cubic in the /c's. 

A particular choice of Cj's define a given Galileon theory. They are considered constant 
over all of spacetime while the values of p{x) that characterize our source will have various 
profiles for different physical configurations. Consider a single point in this space, say Xi. 

In eigenvalue space (what we will call "A;" space) the equations of motion generate a 
surface (or more correctly surfaces — branches — as our equation is a cubic polynomial) which 
depend on p{x). The < inequalities will define volumes in k space; they are independent 
of p{x). Say we restrict ourselves to a particular surface generated by £^ = p/Mp^. 

^It should be noted that "absolute stability" does not imply that all solutions are stable, as there can 
be different branches in our equation of motion surface. But, once again, it does imply that if a solution is 
stable at one spacetime point it is stable over the rest of spacetime. 
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The question is: given particular values of the q's and p(xi) are there some values of the 
fc's that lie on this surface but violate < or Zi < (it is enough to consider Zq and a 
single Zi)7 If this is the case, then either the entire surface is imbedded inside a region where 
at least one > 0, or this surface intersects with the 'marginality surfaces' — the surfaces 
generated by = If this particular surface of solutions fails to intersect with any of 
these marginal surfaces, and has at least one particular choice of fc's where the stability 
inequalities hold, then we say that this surface is a stable surface at the point Xi. 

Repeating this analysis at a different point in space, say X2, means only taking p{xi) — )■ 
p{x2)- Beyond restricting the sources to positive ones, p{x) > 0, a priori we have no idea what 
the source profile, p{x), will be. Thus, in order to ensure stability of a solution generated 
by a particular source configuration p{x), we want to find the {cj}'s such that the family of 
e.o.m. surfaces generated by the possible values of p{x) is fully embedded within the stable 
regions (Z^ < 0). That is, the possible e.o.m. surface at any spatial point is a stable surface. 

Usually, the e.o.m. surface S = p/Mpj will have multiple branches. Consider a particular 
solution TTo that at one point xq in real space is a point (in k space) sitting on a particular 
branch of the e.o.m. surface S = p(xo)/Mpj. If our solution is continuous then at a different 
point Xo + 6x this particular solution will on be the same branch. That is, the branch it is on 
now is obtained from the previous branch under a continuous change of p : p{xo) — )■ p{xo+Sx). 
Solutions are confined to a single branch — they cannot jump from one to another. We can 
therefore analyze each branch in isolation. 

For the purposes of our proof there are two marginal surfaces in k space we will be 
concerned with. The first is defined by Zq = and the second by Zi = 0. We will consider 
their intersections with the S = p/Mpj surface. If it is possible that no intersection occurs, 
then we must check which side (the stable or unstable) the particular surface generated by 
the equations of motion falls. In total, consideration of each intersection will generate a set 
of 'stable' (in the sense given above) {cj}'s for each marginal surface — if they exist. The 
intersection of both these sets will be the stable choices of the {cj}'s for that particular 
S = p/MpY surface. The intersection of these stable choices for every positive p will give the 
absolutely stable values of the {cjj's for the single field Galileon. 

In summary, searching for absolute stability corresponds to looking for values of the q's 
(if any) where some particular branch of any e.o.m surface S = p/M^^ does not inhabit any 
of the volume excluded by the stability inequalities regardless of the value of p. Given that 
the entire e.o.m. surface steers clear of regions of instability, we are assured that if there 
exists a solution whose value at one point in x space happens to be associated with some 
point on the stable branch of the e.o.m. surface in k space then this solution is stable over all 
of space. Calculationally, we determine the absolutely stable region by taking the following 
steps: 1) find the conditions for {q}'s such that there are no real solutions to the algebraic 
equations Z^(A;i, ^2, ^3) = and S{ki, k2, k^) = p/M|j (for a fixed p), 2) check that some 
point on the S{ki, k2, k^) = p/Mpj surface satisfies Z^ < and 3) repeat this process for all 
p e [0, +00) to obtain the intersection. 

^Just for emphasis: the space that all these surfaces live in is, of course, the {ki,k2,k^) space. The only 
function (at this point) of the real space, (cc, y, z), is to give us the single point xi whose purpose is to pick 
out a value of our source, p{xi). 
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3 Details of Analysis 



As previously mentioned, in the presence of non-relativistic sources we can diagonalize the 
matrix dadpiTQ at a point by an appropriate Lorentz transformation. We can then write the 
Z's, Q, and the equation of motion of ttq, ([2]), as 



C2 + 2C3 



^ ] ka j 



va=0 



— 6C4 



1 \ 



+ 24c5 



(14) 



333 3 
£ = ci + C2 ^ /ca + C3 ^ fca/c/j + C4 ^ /e^/^^A;^ + C5 ^ k^kpk^ks - = (15) 

where dadjjT^Q = k^p = diag(fco, fci, ^2, ^3)- As mentioned above, classical stability corre- 
sponds to ^0,^1,^2,^3 < 0. 

Since we have a non-relativistic source we can consistently suppress the ko dependence 
in the above expressions as it is suppressed by f ^ and recover the static limit. In particular 

P 



S = ci + C2{ki + k2 + ks) + 2c3{kik2 + kik^ + A;2/i;3) + Qci{kik2k2, 



Ml, 



(16) 



and 



Zq = - [c2 + 2c3{ki + k2 + ks) + Qci{kik2 + k2k3 + kiks) + 24c5kik2k3\ (17) 

Zi = -[c2 + 2c3{k2 + k3)+6c^{k2k3)] (18) 

and similarly for Z2 and Z3. 

In the following subsections, we will proceed by analyzing the various possible scenarios — 
the S = p/M|[ surface intersecting with the = or Zi = surfaces — independently, and 
then consider the intersection of their constraints. To claim absolute stability we must then 
further take the intersection of the combined constraints for all positive values of p. For the 
moment we will concentrate on cases with a nonvanishing C4, and leave the various special 
cases associated with C4 = for the Appendix. 

An important point that we prove in Appendix [C] if one branch of the equation of motion 
surface intersects a marginal surface, then all others do as well. That is, it is enough to find 
a single intersection for a particular set of q's to rule out absolute stability for that set. 



3.1 Conditions for intersection of the S = p/Mp^ and Zq = surfaces 

(C4 7^ 0) 

Assume for the moment that we are in the stable region of the Zj's, that is Zi < 0, Z2 < 0, 
and Z3 < 0. We can express ki, k2, and ks in terms of the particular values of the Zj's. 
These are 



, C3 ■\/2/i/2/3 , 

ki = - —± -7^ , where 



/, = f{Z,) = -3Z, + 



^-3C2>^-3C2>0 



C4 



C4 



(19) 
(20) 
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The last inequality comes from constraints obtained by analyzing the stability and existence 
of radial solutions of the Galileon theory [1] 



C2 > 

C3 > 
C4 > 
C5 < 



|C2C4 



(21) 



We normalize the field such that it doesn't carry any dimensions, that is [vr] = = 1 
where the brackets mean the usual "dimensions of" and M means "dimensions of mass". As 
the action is dimensionless in Planck units we have that 



21 2-n 



(22) 



In order to compare the free parameters in the Lagrangian we define dimensionless quantities 
from the dimensionful c,'s 



ttl = ci - 



P 
Ml, 



y/cl, a2 = C2, as = 



C3 



"5 = 



_C5_ 

3/2 



(23) 



From (21), it immediately follows that 



«3 > \ -«2 , and as < 



(24) 



Note that while originally we had the whole set of parameters {ci, C2, C3, C4, C5}, a stable 
choice of the these parameters depends only on the choice of the dimensionless parameters 
{ai, a2, as, as}. The C4 dependence disappears because we use it as our units — we measure 
everything in terms of C4. |^ At the classical level, provided it has the correct sign, the overall 
normalization of the Lagrangian does not matter. Noting that a2 = C2 > 0, we can therefore 

work with J 23 I which for brevity of notation we define as {ai, as, a^}. Equivalently, 
we are free to normalize our Lagrangian with the simple choice of C2 = a2 = 1. Either 
way, we are left with three parameters {ai,as,a5}. We are now ready to investigate the 
intersection conditions. 



Inserting our solutions for the fcj's, (19), into Zq and the equation of motion, we have 



Zo = -C2 + 



8a3a5 /I 4a3a5\ 
~9 Is 9~) 



'2^2. 



as 



5t 



y2i 



as 



E = ai — a-i + 



4a 



± 



± 



V2( 



On 



9V2 I 3^2 



u6 = (25) 



uS = (26) 



where 



t = /l + /2 + /3, W=| + | + |, 5 = V/1/2/3 

Ji h j3 



(27) 



^We could have used any other dimensionfuh parameter, but we we find C4 a convenient choice. 
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Using (25) and (26) we can always express u and S in terms of t, which we will treat as 
a free parameter. That is, for some fixed value of t we can solve for u{t) and 6{t) using the 
above constraints. Furthermore, we can solve for /i, /2, and /s via the algebraic equation 



Ft{x) = x'- tx^ + u{t)5{tfx - 5{tf = (x - fi)ix - f2){x - /a) = (28) 
Finally, we can then invert 

f{Zi) = -ZZi + 2al - 3a2 (29) 

to obtain the Zj's. 

'Instability': For a particular choice of the Cj's or (a^'s) such that for some value of 
t, the Zj's are found to be negative implies that a solution could cross into the instability 
region. Thus, there could exist unstable solutions (solutions in the Zq > volume). It is 
not guaranteed that the solution is unstable, but rather that this is a possibility. Using our 
terminology: there will be no absolute stability. Hence the quotes. 

Stability: A stable choice of Cj's or (a^'s) corresponds to there being no intersection of 
the surface of equation of motion and the Zq > volume of instability in the {ki, k2, k^} 
space. Thus, we are seeking a choice of a^'s such that, for any t, the equation Ft{x) = 
cannot have three real roots, all of which must be greater than (203 — 3). This can be written 
in the statement 



{t\t > 3 {2al - 3) , A3[F,] > 0, Ft {2al - 3) < 0, 
and F; (x) > for any x < 2a^ - 3} = 

where A3 [Ft] is the discriminant of the cubic equation Ft{x)=0. 

By virtue of the quadratic nature of Fl{x), the last condition can be further simplified. 
Indeed, the axis of symmetry of the upward-opened parabola Fl{x) = Sx"^ — 2tx + u{t)6'^{t) 
is X = I > 2^3 — 3, so the last condition can be replaced by Fj. {2al — 3) > 0. 



Interestingly, one observes that both roots of solution ( 19 ) yield the same value of u(t)6'^(t) 
and 6^{t) upon which our auxiliary function Ft{x) depends. Therefore, both roots actually 
give the same condition for the stable choice of a's and henceforth we can focus on either 
one. 

In summary, the stability condition coming from demanding that the e.o.m. surface not 
intersect the Zq = marginal plane reads 

{t\t > 3 {2al - 3) , A3 [Ft] > 0, Ft {2al - 3) < 0, and F/ {2al - 3) > 0} = (30) 

Now, as mentioned in our general outline, we still need to check that the surface generated 
by the e.o.m. is on the side of stability {Zq < 0) so that we know we are seeing absolute 
stability as opposed to guaranteed instability. But we hold off on this final check for just a 
moment. 
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3.2 Conditions for intersection of the £ = p/Mp^ and Zi = surfaces 

(C4 7^ 0) 

Assume that we are in the Z2, and Zq stabihty region. Using the expression for Zi given 
by (18) we are free to write the equation of motion as 

£ = c, + C2{h + h) + 2c^{k2h) - hZ, - = (31) 
^8 = ci + C2{h + h) + 203(^:2^3) - = (32) 

Solving the above equation together with the Zi = equation yields the solutions 

-3ai + as±VT 1 
k2,3 = 7 — 7= (33) 

6 — y/C4 



where we have already bothered to normalize everything such that 0:2 = 1 and where F = 
9al + 6 — ISai^s — 3a3 + Saiaf. The plus or minus indicated above means that /c2 must take 
the plus while must take the minus, or vice versa. We pick one. Plugging these solutions 
into the expressions for Z2, Z3 and Zq, all of which are negative given our assumptions, yields 

Zo = - Q^2^ ) (9"i - 603 + 2al + 6^5 - 12010305) < (34) 

(9ai - 9^3 + 4a|)fciycl+ (-3 + 30103 + 03) + v/r(3A;iycI + 03) 

= 2^3 -° ^^'^ 

(9oi — 903 + Aal)kiy/c4 + (—3 + 301O3 + O3) — \/T{3kiy/c4 + 03 



^3 = ^ " ^ ^ ^ — ^ ^ y ' < (36) 



2o| -3 



The latter two conditions imply 

(9oi -9o3 + 4oi^)/civ/cI+ (-3 + 3oi03 + o^) > Vf\{3ki^ + as)] (37) 
A stable choice of parameters corresponds to 



{kl^/c4 I A;ia/c4(9oi — 603 + 203 + 605 — I201O3O5) > 0, r > 0, 

and ((9oi -9o3 + 4o^)A;iv^+ (-3 + 301O3 + 0^)) > yf|(3A;iv/^ + 03)|} = (38) 



which could be further simplified upon taking into consideration some of the information 
coming from ( 30 ) to 

03 - ^o^ - ^(2o^ - 3)^/2 < 01 < 03 - ^al + ^{2al - 3)^/^ (39) 

We will leave the lengthy algebraic analysis that generates the above conditions to Appendix 
[B| to which careful readers are referred. 
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3.3 Stable Region: Local stability to absolute stability 



We must, of course, take the intersection of the stable choices coming from both conditions. 



(30) and (39) 



So far, we have worked out the stabihty conditions for a given choice of a^'s, with the 
external source p being taken as a fixed parameter hke the intrinsic ones (cj's) that define the 
theory. Say that, at the particular spatial point we are working at, the value of p is such that 
the e.o.m. surface £ = p/M|j are completely imbedded in the stability region. As we move 
to a new point in space p will generically change. Thus, we need to do our analysis all over 
again for this new value of the source. The convenience of our method is the details of how 
the source changes are washed out as we consider the whole family of surfaces generated by 
the e.o.m. rather than any particular solution. We lose some information, but we have made 
the problem tractable. We don't have to deal with the functional dependence of our source, 
p{x), boundary conditions, etc. All that matters is the range of values p takes, [pmin, Pmax]- 
The values of p encountered in any kind of astrophysical/cosmological application of the 
Galileon theory will be vast, spanning over 40 orders of magnitude from the average density 
of the universe (~ {lO'^eV)"^) to nuclear density (~ (GeVY). To describe a universe like 
our own in our units we can be free to take pmin — ^ and pmax — ^ oo. 

For a particular value of p, we are given a particular volume in parameter space {01,03,05} 
(in units of C4) by satisfying (30) and (39). Let's call this set of stable parameter choices 
A{p). In order to achieve absolute stability we need to take the intersection of the A(p)'s of 
all possible p's. Thus 

00 

Absolutely stable region of {ci, C3, csj's = |^ A{p) (40) 

p=0 



However, it is easily seen that (39) cannot hold for any source, since as p — ?■ 00, ai 



becomes more and more negative (for a given ci) and eventually fails to fall into the region 



specified by the fixed value of (39). This validates us not checking that the surface 
generated by the e.o.m. lies in the stable region. In summary: From the analysis above, we 
see that for generic sources, no matter how we choose the parameters that define our theory 
(finite values of ci, C3, and C5) the surface generated by the e.o.m. will pierce the marginal 
surface generated by Zi = soiling any hope of absolute stability in our theory. 

Moreover, when examining the various special situations with C4 = discussed in Ap- 
pendix, we find that there is no absolutely stable choice of parameters except for the 
C4 = 0, C5 = case, which is exactly the DGP parameterization. 



4 Analysis of stability regions 

At this point the algebraic relations of ( [30| and (39) don't provide much intuition. As an 
illustration we offer Figure [T} a plot of the stable regions associated with no intersection of 
the Zo = and Zi = marginal surfaces with the e.o.m surface for a particular choice of p. 

The dynamics of the graph as a function of ai are what we are interested in. There are 
obviously several curves that generate the boundaries that we see in Figure [TJ but the one 
that is the most prominent is the vertical line forming the lower left boundary Zi = region. 
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Figure 1: Intersection of stability conditions given by = and Zi = for (from left to 
right) ai = —1, ai = —10, and ai = —100. Red corresponds to stable regions for Zq = 
and blue corresponds to stable regions for Zi = 0. The union of the two regions is the set of 
true stable choices for and a^. 



To get an analytical sense of how the total region of stability depends on our source, we 
investigate the analytical behavior of this line. In particular, we see that it moves from left 
to right as we increase our source (decrease ai). We want to find out exactly how it moves 
as we vary the source. If, for instance, it moves logarithmically then it would not be effective 
at restricting the absolutely stable region (we would have to be more careful when we said 
that p —7- oo) whereas if it goes as some power law we can essentially claim that there is no 
absolutely stable region. 

As shown in Appendix [B] the key constraint in the Zi = conditions is that some function 
(defined in the Appendix) r(ai, a^) > 0. Setting r(ai, 0:3) = we find a particular solution 
asi^ai). Examining this solution, we see that for large values of the source this line moves 
like 



"'^V VJ ^^^^ 

and so we see that for generic sources there are no absolutely stable regions of parameter 
space (in and a^) so that classical solutions of the equations of motion are guaranteed 
stabihty. 



5 Outlook and Conclusions 

After a careful analysis of all the various classes of Galileon theories (different values of 
the Cj's), we find something fairly striking. It seems that only the DGP case is absolutely 
stable — the equations of motion don't intersect the marginality surfaces. Whereas for other, 
more general, Galileon theories the equations of motion do intersect the marginality surfaces. 
This means that the strong statements given in [1] about regions of stability for spherically 
symmetric solutions does not extend to the general case. When generic solutions are taken 
into account the only absolutely stable region shrinks to that of the DGP model (c2 > 0, 
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C4 = 0, and C5 = 0). That is, the classical stability of the non-DGP Galileon theory as 
investigated by P is not as strong as suggested by that analysis. We want to emphasize 
that, while a lack of absolute stability does not guarantee instability (actual solutions could 
lie on the e.o.m. surface safely away from the marginal surfaces — the spherical solutions must 
do exactly this), it does seem to suggest that it could exist. When constructing (non-radial) 
solutions in generic Galillean theory, one needs to explicitly check classical stability via 
some method possibly valid only for those particular solutions. Of course, many physically 
interesting systems exhibit classical instabilities (gravitational collapse for instance). Our 
work suggests classical instabilities could be quite generic in single field Galileon theories. It 
would of great interest to find and investigate these potentially unstable solutions. 
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A Special Case of C4 = 

A.l C5 = 0, C3 7^ 

This is the DGP-like case and has been thoroughly discussed in [16]. For completeness, we 
quickly reproduce the arguments allowing for more general terms. Let's put ko back into our 
expressions, as it allows us to treat Zq and the Zi's in a similar manner for brevity (they 
all take a similar form). If any of the = 0, for arguments sake Zq = 0, while the others 
remain < 0, then the equation of motion can be written 

^ = -^^+Ci-^^-^^{{Z,-Z,r + {Z,-Z.,r + {Z,-Z,f}=0 (42) 



Thus, if 



_c2 + 3c3(ci-^)<0 (43) 



there will be no intersection. In particular, we want (43) to hold for all p > 0. Therefore we 
demand that 



c2 

C3 > 0, and ci < — (44) 

3C3 

Moreover, consider a solution with an known asymptotic behavior, (for instance, at infinity, 
it is deSitter-like: ko = ki = k2 = = k) . For a localized source, it is easy to show that the 
solution 

-C2 + - SciCa 
k as cx) (45) 

DCs 
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is stable, since = — C2 — Qc^k = — a/c^ — Scic^ < 0. 

In the case of a non-localized source, or more complicated asymptotic behavior, we may 



be not able to check its stability directly. However, by constraint (44) we have absolute 
stability and so, provided our solution is stable at one point in space, we guarantee stability 
throughout the whole universe. 

A.2 C3 = 0, C5 7^ 

Zi = —C2 < 0, so the only possibility is that the equation of motion £ = p/Mp^ (for some 
fixed p) may intersect with the Zq = surface. In that case we can rewrite the Zq = and 
£ = equations as 

h + h = -^i^i-l^)-f'i (47) 



We can see that /c2 and are the roots of the equation (treating ki as a free parameter) 



for any C2 > and C5 < (remembering that these are conditions demanded by the analysis 
of absolute stability of spherical solutions) we can always choose a ki such that the above 
expression has two real roots, and thus an intersection between the marginality surface 
Zq = and the surface generated by the equations of motion (for any p) occurs. According 
to our logic, we cannot say such a choice of our parameters leads to an unstable system, but 
rather that there is no absolutely stable choice of q's where C3 = C4 = 0. 

A.3 C3 7^ 0, C5 7^ 

A. 3.1 Conditions for £ = p/Mp^ surface intersecting with Zq = surface 

We can rewrite the Zq = and the £ = p/Mpj equations as 

K.K^Ks = -^(c2 + 2t) (49) 
K^K^ + K^K^ + K^K, = -i(/3i + c2t) (50) 
where we have defined the dimensionless parameters 

K. = csh, 13, = {ci - ^)cs, /35 = ^, t = Ki + K2 + Ks (51) 

Note that 

Z, = - (C2 + 2(7^2 + Ks)) , , Z2 = - (C2 + 2(7^1 + Ks)) , , Z3 = - (c2 + 2(7^1 + K2)) 
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It is easier to discuss a cubic equation (about x) whose three real roots are K1 + K2, Ki + K^, 
and K2 + K3: Gt{x) = 0, where 



Gt{x) = {x - K,- K2){x - K^){x ~ K2- K^) 

= x'- 2tx^ + [{K, + K2){K^ + K^) + {K, + K2){K2 + K^) 

+ {K, + Ks){K2 + K3)]x - [Ki + K2){Ki + K^){K2 + K3) 
= x^ - 2tx^ + [f + {K1K2 + K1K3 + K2K3)]x 

-{K^K2 + K^K^ + K2K^)t + K^K2K^ 

= x'- 2tx' + [t^ - + C2t)]x + + C2t) - ^(C2 + 2t) (52) 

where we used the identities 

{K, + K2){K^ + K3){K2 + K,) = {K,K2 + K^K^ + K2K,){K, + K2 + K^) - K,K2K^ 

and 

(in + K2){Kl + JsTS) + {Kl + K2){K3 + K2) + (Kl + K3){K3 + ir2) 
= (Kl + K2 + K3f + (K1K2 + KIK3 + K2K3) 

If the three real roots of Gt{x) = are all greater than — C2/2 for some value t > — 3c2/4, 
in which case Zi < 0, then the equation of motion surface S = p/Mp^ intersects the Zq = 
surface. Thus, a stable choice of the Pi parameters corresponds to 

t > AslGtix)] > 0, Gt{-j) < 0, and G[{x) > for any x < -|| = (53) 

By the same argument as above, the last condition can be replaced by (?[(— y) > 0. 

A. 3. 2 Conditions for S = p/Mp^ surface intersects with Zi = surface 

Proceeding in a similar manner as above, we can use the same dimensionless parameters and 
the Zi = and £ = p/Mpy equations to write 

K2 + KS = -| 

K2K3 = \{cl-2f3,) (54) 

Note that the Ki dependence drops out of the S = p/Mp^ equation, as the part that is 
proportional to Ki is also proportional to Zi, which by assumption equals zero. Solving 
equations (54), we get 



-C2-A/-3ci + 8/3i; 



-C2 + 



-3ci + 8/3i) 



(55) 
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and thus we can input these values into the expressions for Zq, and Z3. They are simply 

^0 = 
= 



-2i^i(l + 3/35 (c^-2/3i)) 
1 

"2 
1 

"2 



C2 + ^Kx + 

C2 + 4fs:i - 



3c2 + 8/3i 



-3ci 



(56) 



If there exists any such that all these Z's are negative (or vanish), the e.o.m. surface 
may intersect with the Zi = hyperplane. A stable choice of /3j's thus corresponds to 

{Ki\ Zo{Ki) < 0, Z2{K,) < 0, ^3(7^1) < 0, and - 3c^ + 8/3i > 0} = (57) 

Note that Z3 < is sufficient for Z2 < 0. After normalizing C2 to 1, we find that 



If /3i < |, this is a stable region, since it violates the last inequality in (57). 



If I < /3i < |, the Z3 < condition reads Ki > y/S(3i — 3 — 1, which nevertheless 
allows Ki to be either positive or negative, i.e. there exists some Ki such that Zq < 
regardless of the sign of 1 + 3/35(1 — 2/3i). Therefore this is not a stable region. 

If /3i > |, 1 + 3/35(1 - 2/3i) > as /35 < 0. For a sufficiently large Ki > 0, all the 



inequalities in (57) are easily satisfied and so this is not a stable region. 



In summary, in order to avoid the intersection of the e.o.m. surface with the Zi = 
marginal hyperplane, we should demand /3i < 3/8. 

A. 3. 3 Combining the conditions for the £ = p/Mp^ surface not to intersect with 

Zi = and Zq = surfaces 



If we restrict to /3i < 3/8 we^e that the G;(-1/2) = {t + fy + |(| - /3i) > for any t, 



and so the last inequality in (53) gives us no new information. Thus, for S = p/M|j to be a 
stable surface implies 



/3i<3/8, /35<0 



and 



A^lGtix)] > 0, (1 + 2t) [(6/3i - 3)/35 - 1] > 0, and t > - 
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(58) 
(59) 



Now it is time to remind ourselves of the different roles of ci and p, both of which are 
wrapped up in the definition of /3i. Since the requirements of a stable region only constrain 



the intrinsic parameters (the q's), we demand that (58) and (59) hold for any positive p. As 
we will now show, this is not possible for a sufficiently large density (i.e. as /3i — )■ —00). 



Fix /35 7^ 0, for large negative values of /3i, the second and third inequality in (59) gives 



t > Note in particular that t = a/— 2/3i > satisfies this inequality. Now, consider 



A3[Gi(x)] 



;-2/3i 



|5/2 



> 



(60) 
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where ". . . " indicates the subdominant terms for sufficiently negative i.e. we have found 



some t satisfying all the constraints in (59), which implies that for any nonvanishing (3^ and 
a generic source, no choice of parameters {q} is absolutely stable. 

Interestingly, for a vanishing (04 = as well), we recover the DGP-like case, which is 
a stable choice for any positive sources. For a non-zero /^s, stability is only assured provided 
we restrict the source such that 

s{(3,) <(3^< 3/8 (61) 

where ^(/^s) is some complicated function which goes to —00 as /^s — ?■ 0^. It is in this sense 
that the DGP model is special. It is an asymptotic state in comparison to all other possible 
models generated by different points in the parameter space. One has to drive — )■ 0~ in 
order to ensure stability. 



B Further details of analysis 



We will simplify equation (|38|) using some algebraic tricks. In analyzing the conditions for 

(62) 



the Zi = region we begin with the second condition in ( 38 ) 

T = 9al + 6- ISaias - "ial + %axa\ > 
setting r = and solving for ai 



ai = 03 



9 



V2 
9 

V2 



i2al - 3)3/2 ^ ^ 



or «! = as - -ag + —-{2al - 3)^/^ = a+ 
9 9 



(63) 
(64) 



Note, that since 2a| — 3 > both of these solutions are real. Thus, we can rewrite the second 
condition as 

T > ai < a- or «! > a+ (65) 



Now, let's consider the third constraint in (38) 



((9ai - 9^3 + 4:al)ki^ + (-3 + Saias + al)) > Vf\ {3ki^/c^ + 03) I (66) 
^ ((9ai - 9^3 + AaDkiy/c^ + (-3 + 3aia3 + aj)) > and (67) 

((9ai - 9^3 + 4a^)fci + (-3 + 3aia3 + al)) ^ > T{3ki,/^ + 03)^ (68) 



Notice that when ai < a_, the prefactor to ki, (9ai - 9^3 + 4^3) < -V2{2al - 3)^/^ < 0. 
Similarly, when ai > a+ the prefactor > +y/2{2al - 3)^/^ > 0. Thus, we can combine the 



inequalities (67) and (65) as 



ai < a- and kiy/c4 < K{ai,a3) or ai > a+ and kl^/c4 > K{ai,a3) (69) 



where 



(—3 + 3«ia3 + Og) 
(9ai - 9^3 + Aaf) 



"3 



+ 



(2«i - 3)^ 



3 3(9ai - 9^3 + 4af ) 



(70) 
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Now, consider rt68l) which can be rewritten as 



{3-2alf{l-2aia3 + hy/^{-6ai + 2a3) + {kiy/c^f{4:al-6)}>0 ^ (71) 

1 - 2aia3 + ki,/^{-6ai + 2^3) + iki,/^Y{4al - 6) > ^ (72) 

1 3ai — as — vT 1 3ai — ^3 + 
= k2{ai,a3) or ki > 



ki < 



6 



6 



/;;3(ai,a3)(73) 



In the last line, (33) was used. Let's compare (73) with our previous results (69). There are 



two regions we need to concern ourselves with: ai < q;_ and ai > a+. Amazingly, we find 
the combined result is simply 

ai < q;_ and /ci^/q < 1/04^2(^15 0^3) or ai > a+ and kiy/c4> -\/c4/c3(«i, 0:3) (74) 

due to the observation that 

(75) 
(76) 



04/02(0:1,03) < ^/c4k2{a_, a^) = K{a^,as) < K{ai,a3) for ai < a_ 
04^3(01,03) > yc4^3('^+) '^3) = K{a+,a3) > ^^(01,03) for oi > o^ 



In order to have guaranteed stability we need the conditions contained in (74) and the 



first condition of (38), the 05 dependent one, to yield the null set. Equivalently, exphcitly 



separating the two regions, we may write 

• If Oi < o_, for any kiy/c^ < ^04^2(^1, 03), A;iyc4/(oi, 03, 05) < 

• If Oi > o+, for any kiy/c^ > a/c4/i;3(oi, 03), ^1^04/(01, 03, 05) < 

Where we have defined /(oi, 03, 05) = (9 — 120305)01 — 603 + 2o| + 605. Let's investigate 
when these conditions are satisfied - when we are guaranteed classical stability. 

Instability for oi < o_ 



As we are interested in the region where 03 > a/3/2 and 05 < 0, we see that (9 — 
I203O5) > and thus 



/(oi, 03,05) < /(o_, 03,05) 

= (2o^ -3) |-03 - ^JA^^4^+ y(4o3-\/4oi-6 + 8o^-6) 
< for any 05 < and 03 > a/3/2 



(77) 
(78) 

(79) 



So, for oi G (— oo,o_], 7(01,03,05) < and /c2(oi,03) < 0, from which we can see 
that there always exists some ki^/cl such that fciy^/(oi, 03, 05) is positive. Thus, for 
Oi G (— oo,o_], the condition (87) is not a null set and the corresponding parameter 
choice is not a stable one. 

Instability for Oi > 0+ 

Note that for this region ^3(01, 03) is a monotonically increasing function of oi 



^3(0^1,03) > fc3(0+,03) 



-03 + 



03 - 



(80) 



and that A;3(a+,Q;3) < 0. There exists an a^, such that A;3(q;=i,, 03) = 0. In fact, 
a* = 1/203. 

If «+ < ai < a*, then A;3(a+,a3) < A;3(ai,a;3) < ^3(0;*, 03) = and therefore the 
^ia/cI's, satisfying /ciy^ > A/c4fc3(ai, as), could be either positive or negative and so 
be kl^/c4f{al, a^, as), regardless of the sign of /(ai, a^, as). Once again, the condition 



(87) is not a null set and the region ai G [a+, a^,] is not stable. 



If ai > a^,, then /c3(ai,a3) > A;3(a=K,a3) = 0. Since 

(2a2 - 3)2 

/(ai,a3,a5) > /(a*, as, as) = ^ > 

2as 



so it follows that for any fciy^ ( > yc4fc3(ai, as)), /ciyc4/(ai, as, as) is positive. By 
the same argument above, we see that ai > a^, is not a stable region. 

In summary, after a long and torturous process we see that the only stable region (i.e. 



the choice of parameter making condition (38) a null set) is determined by F < 0, or in other 
words 

4 \/2 4 \/2 

a_ = as - -al - ^{2al - 3f/^ < a, < - -al + ^(2a2 - 3)^/2 = a+ (82) 

C Intersection of all branches of equation of motion 
surface S = p/Mp^ with marginal surfaces 

In this paper we have shown that for a fixed p, there are conditions on the parameters 
{ci, C2, cs, C4, C5}, such that the surface generated by equation of motion £ = p/Mpj will 
never intersect the marginal surfaces = 0. However, when p is allowed to vary from to 
00, these conditions generally cannot hold, i.e. no matter how we choose the q's (with the 
exception of the DGP-like case), for some p', the e.o.m. surface will eventually intersect at 
least one of the critical hyperplanes. 

However, since that e.o.m. surface has disconnected branches, it is possible that only 
some branches of the e.o.m. surface intersect a marginal surface, while others are entirely 
buried within the stable region < for all positive p. If there is a chance that this is the 
case, then we cannot claim that any particular choice of parameters is not absolutely stable 
until we investigate each individual branch. 

Here we rule out this possibility by showing that on each branch of the e.o.m. surface 
there are regions where at least one of the Z^'s is positive, for some appropriate positive p 
(with the exception of the DGP-like case). 

Let's first focus on the case with C4 7^ 0. With our usual notation and defining 

Ki = ^4ki (83) 

The e.o.m. surface is given by 

{K„K,,K,)\K, = -6^^^3 + 2a3(ir, + i^s) + l/ ^''^ 
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Generally this surface insists of disconnected branches, due to the fact that for some {K2-, K^), 
the denominator of (84) vanishes. More precisely, suppose (q'l, 5'2, Q's) is some point on the 



e.o.m. surface (for any p), it is in (see Figure [2]) 

Branch 1, if Ggags + 2a3(g2 + gs) + 1 > and > — ^, ^3 > -y 
Branch 2, if Gga^s + 20:3(^2 + ^3) + 1 < 

0:3 tt3 

Branch 3, if 652^3 + 2^3(^2 + ^3) + 1 > and q2 < — ^, gs < -y 



However, if there exists some K2, satisfying simultaneously 

ai + 7^2 + ^3 + 203/^2^3 = 
6K2K3 + 2a^{K2 + 7^3) + 1 = 

some of the branches become connected (see Figure |3|. Indeed, (85) can be solved as 



(85) 



3q;i — 0:3 ± V r 
= 4a|l6 (^«) 

as long as r = 9aJ + (—18^3 + 80:3)0:1 + 6 — 803 > 0. Equivalently 

ai <a3- g«3-^(2«3-3p, or > - -a^ + ^(2«2 _ 3^3/2 ^g^^ 

As p — )■ 00, «! — )■ —00, and so for finite 03 it is easy to check that ^2,3 < — ^3/3, i.e. in this 
limit. Branch 2 and Branch 3 become connected. 
Pick some point {Ki^K2iK^) on Branch 2 with 

= 7^3 = -a,/3, aud lU = ^^f^ ' " (88) 

6^3 — 9 

Evaluating Z^'s on this point, we find that Zi = {2al — 3)/3 > 0; thus Branch 2 cannot be 
wholly contained within the stable region. 

Similarly, consider some point on the e.o.m. surface, with Ki = K2 = = K; in 
the limit p — > 00 or cti — )■ —00, the equation of motion demands K ~ (— It 
is straightforward to show that this point lies in Branch 1. Additionally, at this point 
Zq ^ — 4a5(— ai) > 0, i.e. Branch 1 can not be entirely within the stable region for all 
positive p's. 

Owing to the fact that Branch 2 and 3 will eventually join each other, we conclude (as 
promised) that none of the branches of the e.o.m. surface remain within the stable region 
Zfj^ < for all positive p's. We can also run a similar argument for the C4 = (non-DGP-like) 
case and come to the same conclusion. 
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e.o.m. surface 6=15, with ck3=2 




Figure 3: The surface generated by the equations of motion in 'k' space with 03 = 2 and 
«! = —15. 



22 



References 



[1] A. Nicolis, R. Rattazzi and E. Trincherini, "The Galileon as a local modification of gravity," 
Phys. Rev. D79, 064036 (2009) [arXiv:hep-th/0811.2197]. 

[2] F. P. Silva and K. Koyama, "Self-Accelerating Universe in Galileon Cosmology," Phys. Rev. 
D80, 121301 (2009) [arXiv:hep-th/0909.4538]. 

[3] T. Kobayashi, "Cosmic expansion and growth histories in Galileon scalar-tensor models of 
dark energy," Phys. Rev. D81, 103533 (2010) [arXiv:hep-th/1003.3281]. 

[4] A. De Felice and S. Tsujikawa (2010), "Cosmology of a covariant Galileon field," Phys. Rev. 
Lett. 105, 111301 (2010) [arXiv:hep-th/1007.2700]. 

[5] N. Chow and J. Khoury, "Gahleon Cosmology," Phys. Rev. D80, 024037 (2009) [arXiv:hep- 
th/0905.1325]. 

[6] N. Agarwal, R. Bean, J. Khoury, and M. Trodden, "Title: Cascading Cosmology," Phys. Rev. 
D81, 084020 (2010) [arXiv:hep-th/0912.3798]. 

[7] R. Gannouji and M. Sami, "Galileon gravity and its relevance to late time cosmic acceleration," 
Phys. Rev. D82, 024011 (2010) [arXiv:hep-th/1004.2808]. 

[8] C. Deffayet, G. Esposito-Farese, and A. Vikman, "Covariant Galileon," Phys. Rev. D79, 
084003 (2009) [arXiv:hep-th/0901.1314]. 

[9] C. Deffayet, S. Deser, and G. Esposito-Farese, "Generalized Galileons: All scalar models whose 
curved background extensions maintain second-order field equations and stress tensors," Phys. 
Rev. D80, 064015 (2009) [arXiv:hep-th/ 0906.1967]. 

[10] P. Creminelli, A. Nicolis, and E. Trincherini "Galilean Genesis: an alternative to inflation," 
(2010) [arXiv:hep-th/1007.0027]. 

[11] G. R. Dvali, G. Gabadadze and M. Porrati, "4D gravity on a brane in 5D Minkowski space," 
Phys. Lett. B 485, 208 (2000) [arXiv:hep-th/0005016|. 

[12] M. A. Luty, M. Porrati and R. Rattazzi, "Strong Interactions and Stability in the DGP Model," 
JHEP 0309, 029 (2003) |arXi v:hep-th /0303116|. 

[13] C. de Rham and G. Gabadadze, "Selftuned Massive Spin-2," (2010) 1006.4367. 

[14] J. Khoury, J. Lehners, and B. Ovrut, "Supersymmetric Galileons," (2011) [arXiv:hep- 
th/1103.0003] 

[15] K. Hinterbichler, M. Trodden, and D. Wesleyi, "Multi-field galileons and higher co-dimension 
branes," Phys. Rev. D82, 124018 (2010) [arXiv:hep-th/1008.1305]. 

[16] A. Nicolis and R. Rattazzi, "Classical and Quantum Consistency of the DGP Model," JHEP 
0406, 059 (2004) |arXiv:hep-th/0404 159|. 

[17] S. Endlich, K. Hinterbichler, L. Hui, A. Nicolis, and J. Wang, "Derrick's theorem beyond a 
potential," JHEP 05, 073 (2011) [arXiv:hep-th/1002.4873]. 



23 



